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ABSTRACT 


This thesis deals with computer aided solution to the 
problem of 'placement* of components, and the 'Routing' of wires 
of interconnection between components with . regard to a two 
layer printed circuit board. The main concern of the thesis 
is the implementation (in Pascal) of a force-vector heuristic 
technique for automated placement of components on a PCB and 
integrate it with an existing router with slight modifications, 
as also a concise review of both the placement and routing 
problems with suggestions for future implementations. 

The problem is of importance since miniaturisation, 
in addition to making the wire routing difficult also makes 
wirelength a major contributor to signal degradation when 
coupled with the fact of high switching speeds in circuitry. 
Therefore both the placement problem and routing problem have 
minimisation/reduction of total weighted wire length as their 
objective functions- The routing and placement problems are 
always dealt with separately because of computational complexity 
of the total problem. 

The first chapter introduces the problem and terminology. 
Next two chapters review the placement and routing problem in 
general. The placement problem using heuristic technique, its 
implementation and integration with the router is covered in 
fxirther three chapters. The concluding chapter contains the 
suggestions for future implementation- 
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INTRODUCTION AND TERMINOLOGY 

This thesis deals mainly with Automated placement 
of multiterminal electronic components in a Printed Circuit 
Board Environment, The problem is of importance due to 
two reasons; Firstly the pace of development in integrated 
circuit electronics is rapidly reducing the physical 
sizes of components, at the same time increasing the 
switching speeds available. This makes the wirelength a 
major contributor to signal degradation and the layout 
must minimize the wirelength. Second reason is also due 
to miniaturisation, which makes the wire routing difficult. 
The solution is sought at the placement stage itself, which 
(as an aid to routing) must also satisfy the connectability 
req.ulrements at the same time that it tries to minimize the 
wirelength. 

A two layer PCB consists basically of a ‘Board*, i.e. 
the area allowed for placement of electronic components. 

Each component, whatever be its internal mechanism, provides 
some external pins. Depending upon the circuit function 
different subsets of these pins are to be inter connected. 
The wires of such connection contribute to the total 
wirelength of the PCB. The wirelengths may be weighted. 



!Ehe production of the PCB involves making an 
artmaster which in tiurn is to he made from a scale drawing 
of the components and the electrical inter connections. 

This project is an attempt to make such a drawing with 
the aid of computer, given the description of the hoard, 
the components and the inter connections with 1he objective 
of minimizing the total wirelength. 

>Jhether such a drawing is attempted with the help 
of the computer or manually the basic problems confronting 
the designer are the same, namely ’PIACEMENT’ AND 'R0UTIN&’, 
The two are inter related but treated separately because 
of the inherent compu’"8.tiona-. complexity of the total 
problem, 

1 , 1 TERMINOLOGY 

'B oard * ; A board is the physical area for placement 
of components of the circuit and for accommodating the 
wires inter connecting them. It is defined by its span in 
A and Y directions or by a set of corner points. 

' Modules * : A module refers to aiqy element that needs 
to be assigned to the board with a specific location, i.e. 
a specific X and Y coordinates. A module may have several 
pins that are electrically isolated but bound physically 
together by the module. All pins common to a single module 
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are -treated as the same pin. As a component, the term 
module may refer to any of these: Package, Chip, IC 
or Connector. 

' Nets’ or * Signalsets * : Ihese describe the inter 
connections of modules. A Net consists of all electrically 
common pins on a board. Each net may employ its own wiring 
or inter connection rule specified. 

1.2 CH/iPTER¥ISl ORGMIZATION 

The following five chapters deal with the subject 
matter as follows; 

Ch apter 2 : .An overview of the placement problem, Ihe 
placement problem is defined and three approaches to solu+ion 
are discussed in essence. 

Chapter 3 : An overview of the routing problem. Brief 
discussion about which of the pins to be connected in what 
order and how. The emphasis is on Lee’s algorithm for 
actual routing between pin pa,irs. 

Chapter 4 ; Discussion in detail of a force-vector 
heuristic technique of placement. The implementation (in 
Pascal) of this technique is the main concern of this thesis. 

Chapter 5 : The pascal implementation of the placement 
program is discussed from user's point of view and documented 



4 


Chapter 6 ; An existing routing program [ 10] is 
documented for ease of user. 

Chapter 7 : Conclusion. 



CHAP 2ER 2 

THE PIACEMEHT PROBLEM OVERVIEW 

2.1 PROBLEM DEEIJIETION 

Eormally, the niodu3.e placement problem consists 
of finding the optimal placement of modules in the board 
with respect to some norm defined on the interconnections. 

§|ince interconnecting wirelengths is a major 
contributor to signal degradation, especially in an 
environment of miniature components, it must be reduced. 
There are also a n\imber of other circuit goals: minimum 
wirebuildup in the routing channels, elimination of signal 
echoes, minimum heat dissipation etc. Since for practical 
purpotses it is impossible to directly incorporate all of 
the circuit goals, the total weighted wirelength is used as 
an approximation in trying to satisfy them enmasse. 

Apart from reduction of wirelength the placement 
should also have as its another primary aim as that of 
aiding the wire routing. Therefore in general we can state 
that the placement problem is that of assignment of M modul 
to N sites that satisfies the connect ability requirements 
and also minimizes the total interconnecting distances. 

2.2 APPROACHES TO SOLUTION 


Because of the complexity and magnitude of most 
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practical problems no guaranteed method for optimum solution 
exists. Breuer's formulation of the problem as an integer 
linear programming task [2] provides an exact solution but 
the constraints imposed make the solution of less practical 
value. Therefore methods based on heuristic rationales 
are employed.. They may be classified as follows. 

’ Constructive Techniques *; Useful for initial placement. 

A placement configuration is produced only upon termination 
of algorithm.. 

’ Iterative Techniques ' ; Useful for placement improvement. 
At every stage there is a complete placement a.vaila.ble. The 
algorithm improves it by repeated modification. 

' Branch and Bound Tech n iques ': These are also 
constructive techniques in essense. 

2.3 CONSTRUCTIVE METHODS 

In this placement configuration is formed by adjoining- 
modules to a subset of already placed modules- The input 
data consists of details of borrd, nets, module.? and other 
control parameters and iteratively selects one of the 
unplaced modules for inclusion in the set of placed modules. 
Once placed the modules are not moved. 

These methods have the advantage of requiring 
relati-vely small amount of computer time and yet sufficient 



7 


for many applications, 'Uie specific methods under this 
category are defined by the particular rules of ’selection' 
(i.e, ordering of modules) and 'positicnin,g ' of modules. 

The known methods are; 'Pair~lihklng Methods' and ’Cluster 
Development Methods' due to J.M. Kurt zb erg. 

2.3.1 SELECTION RULES 

To determine the order of module selection some 
measure of connectivity of the modules is required i.e., 
some notion of how 'strongly' unplaced modules are bound to 
already placed modules, A natural measure of connectivity 
between two distinct modules is sum of weights of all nets 
in which both these distinct modules are members. Let us 
call that measure for two distinct modules K and Y as 
P(.i,Y). 

One selection 'function' may be based as follows. 

It will, consider each module L included in the set of unplaced 
modules and each module Y included in Hhe set of placed 
modules and select that module K for placement which has 

P (K,Y) . A second rule may incorporate the connectivity 
from an unplaced module li to all placed modules Y and select' 
such which gives max^P (X .Y ) . 

Yet another selection rule may be based on the 
number of nets to which unplaced module X belongs and select 
such X which has maximum of its nets in the process of 
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placement. A net is considered in tie process of placement 
if atleast one of its modules but net all are already placed, 
A fourth ru3.e may improve upon it and select such X whose 
nets together contribute the maicimum nximber of already 
placed IIlOd^lles. Yet another rule may improve it by taking 
the number of modules in each net of X into consideration 
in the same rule. 

In general a selection rule does not yield a unique 
module. In that case a 'tie-breaking' rule is necessary 
which can be; (a) look-ahead small number of steps; 

(b) Optimise on a secondary function or 

(c) Choose arbitrarily. 

The above :rules define the generic types but many 
other meaningful rules can be defined. The particular 
application, interdependence of net and computer time etc. 
dictate the choice of a selection rule. 

The selection rule i described are all independent 
cf the positioning rule. Hence the modules can be completely 
ordered, if desired, before the positioning rule is applied. 
But when the selection rule depends upon the relative posi- 
tioning of the modules on board the computation becomes 
more complex. 

Note that when a module has all the modules of all 
its nets unplaced, that module is not a candidate for the 
next iteration of selection :nile. 
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2.3.2 POSITIOmFG RULES 

In positioning a modixle, we try to find the best 
slot for it relative to the nodules already positioned. 
Instead of examining all empty slots, only a comparatively 
small subset of slots (’candidate slots') are considered. 
These subsets may be accumulated along with the positioning 
of each new module or an appropriate set may be regenerated 
for each positioning. Eor example, a. simple pattern of 
four adjacent slots of each newly placed module can be added. 
As each slot is occupied it is deleted from the candidate 
set and the slots already occupied are never re-employed 
and unusable slots are never added to the set. In this 
pattern the periphery of positioned modules is always the 
candidate set, H. Ereitag and M. Hanan of IBM have used 
a variation of this technique. 

With each candidate slot accommodating the modu3.e 
to be positioned, we compute the wireleng+h (cost) of all 
the (partial) nets of which this module is a member, based 
on that particular nst^ inter connect icn rule. The module 
can then be placed in that caididste s3.ot which gives tie 
minimum ccst. 

In practice this can no-; always be done due to 
two reasons; one is excessive computer time and the other 
is that since not all modules i.n each of these nets have 
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been placed, the optimal interconnection can not be 
determined at this point. Thus an approximation rule is 
frequently used. Por example, if the rule is minimum 
spanning tree then for each candidate slot we may compute 
the distance to the closest raodule in each of the partiall.y 
placed nets considered and we choose tha,t tslot which gives 
the minimum total distance. 

2.4 ITEHATIi/E METHODS 

This project concerns the implementation of a heuristic 
technique of applying vector forces to component modules 
and iteratively relocating them to arrive at a hopefully 
optimal placement. '2his iterative method, implemented in 
rascal, is based on a IEI;,A.GH Design Automation Workshop 
paper [2]. This method which is independent of initial 
placement will be discussed in detail in Chapter 4. Other 
iterative methods which have been tried mainly as an 
improvement of initial placement are: 

(a) Steinberg's Assignment Method 
(Steinberg and Rutnam 1964) 

(b) Relaxation Method (Casey and West 1967) 

(c) Pairwise interchange Method (Kurtzberg and Hanan 
1970) 

(d) Stochastic Method or 'Montecarlo' Approaches 
(Cooper 1964, Kurtzberg I960). 
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2.4.1 iffi STRICT ALGORITHM 

This class of algorithms in general , abstract sense 
can be formulated in the fol3.ovTing paragraphs. 

Step 1 : Generate initial placement configuration. 

Step 2 °. Generate a 'transformation' whereby the 
la,yout is mapped to a new configuration. 

Step 3 i Compute the cost of this new configuration 
and compare it with the old configuration. 

Step 4 ; If the cost (may be wirelength) of the new 
layout is reduced the new configuration is substituted 
for the old one. Otherwise the old placement is retained. 

Step 5 : Apply 'stopping rule' with outcomes j 

(a) Go to step 2. i.e. continue modifying the placement. 
The basis to terminate this loop depends upon 

the rate of improvement in the cost or upon the 
predetermined number of transfozmiations without 
improvement. 

(b) Go to step 1. i.e. a new starting poiht is made 
with a new initial placement for modification. Hie 
decision to stop this loop depends on the 

time available. 


(c) Stop. 
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Ihe resultant placement for each initial placement 
is 'locally optimum' in the sense that there does not exist 
a placement with a smaller cost for the given initial 
configuration. 

We can generate a large number of locally optimum 
placements, thereby increasing the probability that atleast 
one of these placements is close to the global optimxjm. 

Alternatively we can find better locally optimum 
placements thereby increasing the probability that each 
initial placement leads to a placement close to the 
optimum. 

2.5 BRMCH-AND-BOUND METHODS 

011636 methods have the feature that they can be used 
to find an optimum solution (when the placement problem 
is viewed as a quadratic assignment problem) if the order 
is less than 15 modules. Gilmore (1962) and. Lawler (1965) 
have independently applied ''liese methods with Gilmore 
modelling an approximation to the exact branch and bound method. 

Briefly stated, the branch-and-bound method proceeds 
as folO.ows. The set of all feasible solutions is partitioned 
and a search for the optimum is made in each partition. 

1 lower bound is computed for the solutions in each partition, 
and the search in any partition is terminated upon the lower 
bound exceeding the co.st of some previously found feasible 


solution 
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The. feas^ihle solution may be obtained by any method, 
e.g., by a constructive method cr by b ranch- and-boimd 
algorithm itself. 

To define methods for computing lower bounds and 
partitioning sets of lower bounds consider the quadratic 
assignment problem in which there is a connection matrix 
C = and a distance matrix D = [d. .] and we seek a 

permutation p such that 

'^P(i) ^(j) 

is a minimum. 

As a preliminary consider the branch end bound, method 
from the view point of an (n~l) stage decision process. At 
each stage k, a module i is chosen from the remaining (n-k) 
modules and assigned to all the remaining (n-k) slots. The 
term branch and bound comes from viewing the method as a 
search in a decision tree. 


On the first level tree jn Fig. 2 , 1 , module i has been 
chosen and each vertex at first level is a possible slot 
for this module. On the sec.'iid level module ±2 is chosen 
and assigned to any of the remaining slots. 


For each vertex in the decision tree, a lower bound 
is computed. Gilmore developed two lower bounds for 


c . 


ij-' P(i)p(j) 


based on the observation that if 
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c = (c ,Cp,...c ) and d = (d.rdp,...d ) are two vectors 
then the minirnmn dot product of c and d i.e. ISh-c^d^^^^ 
over all permutations p is obtained by arranging c in 
increasing order and d in decreasing order, and then 
multiplying term by term. 

A simple lower bound for the quadratic assignment 
problem is obtained by treating the nondiagonal elements 
of matrices c and d as vectors of length n(n-1 ) and 
performing the proper ordering and multiplication. Diagonal 
elements are eliminated since no pin is connected to itself 
and txiTO modules can not occupy the same slot. Also if some 
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modules are already placed, then the contribution of these 
placed modules to the total distance is computed, the 
corresponding entries from c and d matrices are deleted 
and a new lower bound is calciolated by forming vectors 
from the matrices without the entries used in the computation 
of placed modules contribution. 


In another calculation of lower bound the problem 
reduces to finding the solution for assignment problem 


of finding the minim-urn of over all permutations p 

vrhere [a. .] is the ’assignment matrix’ A. The element a- ^ 
of matrix A is the minimum do'j product of vectors and 
■which are formed from matrices c and d after deleting entry 
'^ii ^3 3 respectively. Since the total cost of placing 
module i into slot 3 i s lower bound 

on the cost of placing module 1 in slot 3. A permutation % 
for which a.„/ . \ is minimum then yields a 3 .ower bo\md for 

XU ) 

the original quadratic assignment problem. 


The efficiency of the branch and bound method depends 
on both the bounding and branching techniques. The function 
of the lower bound is to prune the decision tree; the better 
the bound the sharper the pruning. However a tradeoff is 
required since better the boimd, the more 
computation time required to compute this bound. 


the 
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ROUTING PROBLEM 07ER?IEW 

The process of specifying the precise conductor 
paths necessary to properly interconnect the module pins 
of each net subject to the imposed constraints of wire 
thickness, spacing, constraints, feedthrough reduction 
etc. is the crux of the problem. Our placement techniques are 
aids to this,. We may say that the constraints stated »bove 
contribute to the weight of wirelength which we seek to 
reduce. Again a compromise between running time and optimality 
of solution must be effectedi Various authors agree that 
typicalJy, a ’poor’ solution could be obtained in about 
5 seconds of computer time and the ’best' solution in 
probably over 5 hours. The basic approach to the routing 
problem is stated in the following paragraph. 

'lee’s algorithm' is the basic tool in the final lay- 
out process. But a successful layout results from a series 
of closely interlocked steps, each designed towards the 
ultimate goal of making the final laj^'Out procedure as simple 
and as efficient as possible. The four basic steps are: 

(a) Wire list determination; 

(b) layering; 

(c) Ordering; 

(d) Wire layout. 
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Rougjily speaking these are the *what ' , where', 'when* and 
*how' of the problem. 

5.1 WISE LIST DETERItINATIOI 

This lists precisely what wires are to be connected, 
i,e. which pin to which pin. This is simp3.ified if such a 
list is produced in the output of placement algorithm itself 
or it can be incorporated during the routing process. The 
desired wiring rule for each net, •when Implemented keeping 
the objective function inview produces a list of such pins 
or pin pairs for each net to bo connected. 

5.2 LAYERING 

This involves tentatively assigning each wire to one 
of the layers of ‘the given multilayer board. This is easily 
and arbitrarily solved in the case of two layer bo^^rds. 

One approach is to try to route as many wires as possible 
on the first side under consideration. Another approach is 
discussed in [s]. ■ 

5.5 ORDERING 

This is to determine in xfhat order should the wires 
on each 3.ayer be processed. VJe may consider each wire to be 
a net in the sense we defined a net. Ihen it becomes a 
quest! n of in what order each net is going to be laidout. 
There are many rules of thumb developed for this and there is 
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nothing sacred about a parti cuJ.ar rule and most of them 
tend to lead to a nesting of vfires. Ihe significant 
point is that of finding paths virhich will cause the least 
diff icu3.ties for subsequent path layouts. Ihere is no 
firm guarantee or precise answer to this. However, a 
reasonable choice seems to be that path which uses up the 
fewest grid segments, 

3.. 4 WISE ljQrOUTSl.EB’S ilGORIfHK 

The previous three topics have been specifically 
oriented towards how the individual net wires are to be 
laidout interconnecting each net, Lee's algorithm is a 
funaamental tool which greatly simplifies the process. 

Oliis is actually an application of the shortest path 
algorithm used in O.R, and graph ilieory. Here the paths 
under consideration must foll-ow" the cells of a rectangular 
grid and hence the algorithm becomes' even simpler to apply. 

The process begins with the selection of either pin as 
a starting point. Say A is selected{with reference to Fig.3.l). 
Initially a 1 is entered in each empty cell immediately 
adjacent to the cell containing pin A, Next 2s are entered 
in all empty cells adjacent to those ce3.1s containing 1 ’ s. 

Next 3’e are entered adjacent to 2s and so on. 

The process continues in this manner until one of two 
results occur: Either all paths become blocked, i.e. on the 




Pig. 3.U Lee's Algorithm 

kth step there are no empty cells adjacent to those 
containing (k-1 ) or the target cell B is reached. In this 
case B is reached on the 13th step indicating that there 
is a path from A to B and the shortest such path is of 
length 13. All. that remains is to actually label this path, 

3.4.1 BACK TRACKING 

.Since B was reached on the l3th step, it fol'’.ows 
that there must he a 12 immediate3.y adjacent to it and which 
inturn must he adjacent to i1 and so on. Proceeding in 
this manner we lahel the paths leading hack to A. Now 

' f" -1 ' 

when we have a choice in hack tracking, i.e. when there are 
■■two or more cells with (k-1 ) adjacent to a cell containing k. 
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in theory any of these cells may he chosen and path still 
foimd. In practice (reduce bends) we adopt the guideline; 
Don’t change directions unless you have to. A3.so when 
a txirn must he made we select and try the same sequence of 
directions for all turns to try to achieve a uniform 
nesting. More sophisticated nesting techniques have been 
tried but the program gets more and more complex. 

5.5 SOME .SPEEDUP TECHIIQUES 

In a typical wire layout, 70 to 80 percent of running 
time is spent on Lee’s algorithm and any technique that can 
be employed to speed up its execution has observable pay 
offs. Three commonly used techniques are: 

Starting point selections : Consider the pins I and Y 
in Fig. 5.2. If X is chosen as starting point (Idg. 3.2(a)), 
less than half the total cells will have to be labelled 
before Y is reached. On the other.:’ hand choosing Y as ^ 
starting poini leads to labelling practically all the cells. 




Fig. 3.2(a) Fig. 3.2(b) 

Fig, 3.2: Starting Point Selection 
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So the guide line can hej start from the pin farthest 
from the centre of the hoard. But this loses importance 
as the number of wires laid out increases. 

Double f an out ; In this we begin at both pins. Ihe 



Pig. 3.3; Double Fan out 

program will be much more complicated. 

Framing : in artificial rectangular boundary is 
imposed about the pin pair being processed and no labelling' 



Fig. 3^4; Framing 

is allowed beyond this boundary, say about l0 to 20 percent 
larger than the rectangle defined by the pins. \'Jh.en 
path is not found in the frame it is repeated with 
larger frame. If the repetition occurs too frequently 
towards the end the framing is discarded, 





22 


5.6 ST0H_4,.G-S PROBLEMS 

The primary disadvantage of lee's algorithm 
is that fo'' eveiy cell in the grid there must he a 
corresponding memory location. In certain modified 
algorithm it has been tried to avoid by having a 'list 
of current latest cells' and in each new step the list 
is updated xfith only those cells adjacent to cells on 
this list. Another economy is by economising the label 
each cell has to receive. ’''Oien path length is long the 
label number is more, and more ?.lsts are required. This 
can be avoided by using a b;lnsry coding of only 2 bits 
of 3 bits for each label by utilising the fact that during 
the back track process, it is only necessary that 'predecessor' 
labels be distinguished from 'successor' labels. 
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A FORCE VECTOR HEURISTIC TECHNIQUE FOR PLACEMENT 

This plo.cement technique has been implemented 
in pascal* A listing of which is attached as Appendix II, 
This is a. heuristic technique of applying vector forces 
to modules and iteratively relocating them. The methods 
used to compute these forces improve the ability to find 
solutions that are independent of initial conditions and 
near optimal. This employs a. very flexible modelling 
3 cheme . 

The aims of this technique ares 

(a) Flexibility in handling many sizes and shapes 
of components. 

(b) Adherence to any wiring (interconnection) rules 
specified. 

(c) Minimization of total wirelength. 

(d) Improving connectability. 

(e) Handling practical problems without requiring 
impra.ctical memory, 

(f) Fast and computationally efficient method. 

4,1 MODELLING CONSIDERATIONS 


The ideas of Board, Modul.o and Net are retained in -the 
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same sense as already defined. But there are six classifi- 
cations of modules and three interconnection rules 
implemented. 

4.1.1 MODULE TTPES 

Free modules : 'Siose modules that are free to assume 
any X and Y coordinate necessary for their placement. Only 
the boundaries of the board and interference of other modules 
restrict their positioning. ’.'OIig majority of the modiiles 
are usual3.y of thj.s free type. 

Fixed modifies : These modules can not be moved. Their 
X and y coordinates are specified in ihe input data and 
rema.in so. The terminal connectors for example are of 
this type , 

Fixed~X modules ; These are modules whose X~ coordinate 
must remain constant but may assume any specific Y coordinate 
during the, execution of the pi’Ogram. 

Fixed-Y modules : Similar to fixed“X, But here the 
y coordinate must remain fixed. 'The program may alter 
its X~coordinate . Eorraally seme stand-alone discrete 
component or other special circuit requirement may give 
o ccassion for a fixed X or fixed Y module. 

X-bar modules : The X-bar modules are like the fixed-X 
modules only as far as the fact ths,t their X coordinate must 
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remain constant. But their Y coordinate instead of being 
a single point appears as a vertical line of indefinite 
length to the program. Any connection to such a module 
is niere3,y a perpendicular to that line. 

Y-bar module s ; These are idle counter parts of Z-"bar' 
modules with their Y coordinate constant and appearing 
as a horizontal line of indefinite length to the program. 

Bar modules generally apply to a set of identical components 
occupying a certain fixed area and the connections to which 
f .-om any point of the board is likely to be a. perpendicular 
line to that general area. 

Every module has a physical size. Sizes are 
rectangular dimensions of length and width from the centre 
point i.e, coordinate position of the module. For the 
■placement program modules are modelled as single dimensionless 
points (i.e, coordinate position) and pins of a. module 
are al.so treated as the same point. 

The terminal connector plug can be 'segmented' to 
avoid too manjr modules that would arise by treating each 
board terminal pin as a module. 'Segmenting' is simply 
grouping pins of the connecter plug into areas that are . 
approximately the size of the other modules oh the board. 
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To model irregular shape of bosrds extra sites 
may be defined by adding fixed ts'-pe dummy modules, The 
dmmy modules have no connections and no assignment would 
be made to them, but will exert a physical force on t: e 
positioning. 

4.1.2 BXMPLE OP BOARD MODELLING 

In this board the origin of coordinate system can 
be assumed at left most bottom corner. Each integer 



Pig. 4.1: Circuit board to be modelled 

coordinate unit along the axes may correspond to 
ore 1 or 2 mm as desired. The coordinate positions 
(centres) and sizes of all modules are in integers. 

Modules 1 to 1 1 are of free type with other associated 
information. Module 12 is a fixed type dummy module with 
no connections. Modules 13 to 17 are identical and can be 
modelled as a single X-bar module with no pins though its 
sub-components have pins. The plug connector (pins 19--57) 
can be segmented as three fixed modules with 6 pins each. 
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4.2 SET IN'TSRCCNKEC^ION RULES 

Three basic wiring rules are briefly mentioned. 

Star connection ; This connects all the pins in the 
net to the source (specified modii2.e-pin) directly. 

Serial connection : Connects all the pins in a net 
in a serial fashion by creating a single string from the 
source pin to the final sink in an optimum manner. 

Minimum distance : Connects all pins in the net in 
an optimum fashion without regard to source or sink or 
fan~out restrictions on a pin. 

In addition to these above wiring rules the designer 
may specify any unique wiring rule in a free standing 
subroutine. Critical wirelengtiis can be monitored. 
Artificial nets may be added. The module are also modelled 
very flexibly. Therefore by using any clever scheme 
almost ah 7 design problem can be modelled. 

4.3 PROGRAM OVERVIEW 

This is explained in the flow chart of Pig. 4.2. 

The input data ; This is data necessary to describe 
the model to the program. Consists of board dimension, 
number of nets, .modules, wiring rules, module descriptions 
etc. and other program control parameter. 
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Fig. 4.2: Flow chart of Force Vector Method 


«f! on verge* process ; An iterative process that 
incrementally moves the modules to an oblique layout. 

Modules are considered dimensionless points within a 
system of forces. Also reorders net interconnections 
between iterations in an attempt tc optimize them based 
on wiring rule . 
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'‘Expand* process : This has not been implemented 
in the pascal program of this project. This process 
considers modules with true physical dimensions. Physical 
overlap creates a new type of force that acts to separate 
the components. Electrical attraction forces still app'^y. 
Boundaries of board contribute additional forces. 

Final *assignment * ; DTeceosary for the assignment of 
modules that must be mounted on predetermined sites, i.e. 
orderly, rows and columns. Becoming increasingly important 
with use of iCs and automated production uses Hungatrian 
assignment algoritian. 

4.4 COlfFERGE PROCESS 

The aim of this process is two fold: An optimum 
orientation for all the modviles with respect to each other 
and an optimum interconnection for each net. Optimality 
is with respect to the objective functi n of total line 
length. This iterative procedure has two steps in each 
iteration, which infact is a coordi.nate transformation of 
modules: Firstly all forces are computed for all modules 
in their current position, i.e,, in the same state of the 
process. Secondly all modules are then relocated by a 
fractional amount of the resultant vector at the same time. 
The incremental approach eliminates dependSv^c^ upon the 
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order of selection of modules for the move. The incre- 
menta3. fraction (i.e. ’Mobility Pactor') is chosen to 
give speed of computation but at the same time it should 
not destroy the general orientation which was the basis 
to ca3..culate this move, 

4,4,1 TSCTOH POROES ON A MODULE AND SORT-DELETE 
TECHNIQUE OP COMPUTINl .RESULTANT 

The process begins by computing the best inter- 
connections for each net according to the desired wiring 
rule and the current module layout. 

Each line (i.e., the rectilinear or 'Manhattan' 
distance) between a module and its connected neighbour 
is treated as a separate S vector and Y vector on the 
module. 

Starting with the vector pairs of greatest 
magnitude (X and Y separately), pairs of vectors are 
deleted from the module. All rejaaining vectors (X and Y 
•separately) on rhe modu3.e mast be of same sign and 
they are added together into one vector. Figures 4.3 and 
4.4 illustrate how this technique (named 'Sort-Delete 
Technique') operates and how it minimizes the rectilinear 
distance (which is what we seek) rather than the 
pythogorean distance. 
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Fig. 4.3(a): Module with all 2 vectors equal to linear 
2 distance to connected modules.. 



Fig. 4.3(t>): One pair 'f vectors of opposite sign 
(greatest magnitxide) deleted. 



Fig. 4.3(c): Re suit ant ve c t o .r 

Fi g. 4.3 Sor t-Jelet e Te chn ique 



Fig. 4.4(a) Module X is in equilibrium for simple vector 
sum. 

O (5^^ 

''O 

Fig. 4.4(b); Resxiltant vector in sort-delete technique 
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o 

Pig. 4.4(c)s Module X stable- in Sort-delete technique «• 
Rectilinear distance 

Pig. 4,4 Reduction of Rectilinear distance in 
Sort-delete technique 

4.4.2 KERUI-SIOR 1/EC TORS 

These forces are also to be considered scttng on 
a module before computing i±ie resultant vector on a module. 
This vector arises between the unconnected modules. 

These repulsion forces prevent the electrica.1 attractions 
from pulling all the modules together in a central mass 
and also improve the stopping point by contributing to ensure 
that modules are on the proper side of each other at the end 
of converge process. This is iJJ.ustrated in Pig. 4.5. 

The amount of repulsion is an input variable. The 
paper suggests distance equal to two times typical module 
sizes. I'/hiie this is retained in the implementation 
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it was found useful to drop the repulsion forces 
altogether once unconnecoed modules are separated 
by the same distance (twice typical module size). 



Fig. 4.5(b) Overlap between unconnected B and C 
at this stopping point. 



Fig. 4.5(c) Repulsions between B and C but still their 

length to A is the same. 



Fig. 4.5(d) Better stopping J;paj 3 .ditions with same minimum 
line length 
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4.4.3 MODULE TYPES AMD YSCTOHS 

The vector calculations must take into account the 
differences in modelling tyx^es. Vectors in X-directton 
•^nd vectors in Y-direction must be treated separately for 
each combination of module types- Guide lines are; 

(a) Free modules ; Normal calculations of vectors 
when the other modiile is a free module. If the 
other module is a fixed module then the free 
module receives tvrice the vector force. 

(b) Fixed modules : Receive no vector force. 

^-bar and Y~b ar modules : Treated as fixed 
modules- In the case of X~bar, the other reacting 
module will receive twice the force only in 
^-direction and no force in the Y-direction. The 
case is reverse for Y-bar, 

(^) 1 " Sliced and Y , - FiX^d jaodul Treated as fixed modules 
only as far s.s l~vector 3 and Y-vectors are 
concerned respectively. 

4.4.4 MOVEMENT OP MODULES 

Once the resultant vectors have been computed for 
all modules, the program adds the resultant vector times 
some scalar coefficient to the -position of each module 
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which is its new po^itj.on. Jliis scalar coefficient is 
cali.ed the 'MOBILITY BAG TOR' and is calcnl-ated automatiC8JJ.y 
according to the formula; 


Mobility Factor 


-1 (Average Vector) 
2 ( M 8 X iniiS V e c t or J 


The one-half in the formula is due to the fact that 
vector computed between aiodu3.es J and K. is applied to both 
modu3.e J and module K and if' no other module were included 
in the calculation, modules J and K will meet exactly half 
wa.y between them in one iteration . 

Secondly the formu3.a should help in driving the 
solution to a stable finish for the following reason. 

At the beginning of the carsrerge process, the placement 
may be very poor end vectors computed may be very large. 

If '"he maximum vector is large, the mobility factor 
will be smaller thus providing a desirable initial control. 
As the process continues s.;.id the forces decrease the 
mobi3.i:-y factor also increfses as the vectors approach 
average value. This increase of mobility factor 
discourages the process 'from becoming slower as more 
iterations are completed. 
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4.4.5 HET INTBRCOMECTIOH RUXES AFD ITS FREQUENCY 
OP IMPLEMENTATION 

There is no way of determining the optimum order 
of interconnection for a particular net without considering 
all other connections on the board. In fact until a final 
placement is made there is no way of determining what effects 
all the modules and all the nets will have on a particular 
net interconnection. The wiring rules we implement for 
optimizaticn apply only to any single state of this system 
of forces and moving modules. 

In converge process, we can reorder (i.e. implement 
the wiring rules) the net between any +wo iterati'^ns. This 
however would be time consuming. It ivS apparent th'^t 
between iterations involving very small changes in line 
length, no reordering woul.d be necessary. However, when 
line length is changing sha.rply from one iteration to ihe 
next the reordering would bo required. The program should 
therefore monitor this line length (which it does for 
determining the stopping conditions also as we shall soon 
discuss) and choose pan’ticular frequency of iterations 
after which it will reorder again. This comparison can 
be made on the first iteration after a reordering of net. 

This reordering frequency (in terms of iterations) is a 
matter of experiment and experience. The paper suggests 
the following function depicted in Fig.. 4.6. 



4.4.6 PROGIffiSSIVE ADDITION OF NETS MD INTIAl CONDITIONING 


Why progregsive addition ? The numerous combina'tions 
possible tend to add to the already numorcus possible 
states that could be stable. Therefore, to help insure 
that the process approaches its stable state in ihe most 
unbiased increments possibj.e nets are included progress- 
ive3.y in the process. 
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Which nets first? A characteristic of mos'^ circuit 
designs is that most nets have just a pair of pins, or onl^y 
three modules (pins) per net. Because there is only one 
way to connect pin pairs and that they also foim the 
majority, it is essential that these nets be organized 
first. CONVERGE, therefore, begins by moving all modules 
to a stable state by considering only two pin nets. Once 
a stable state for forces generated with these nets is 
reached all nets with three or less pins are included and 
so on imtil all nets are included. 

I n i t i al. cond iti onin g ; In an effort to further 
improve the \mbiased starting condition of the problem 
all runs begin with the ’star' rule of Interconnection 
regardless of their wiring rule to be implemented. But 
this is only for the First Pass through CONVERGE. 

4.4.7 STOPPING CONDITION 

This is a point of relative stability where total 
line length stops changing to within 'some small value' 
between iterations. This value again is a matter of 

■ -7 

experiment and experience. The paper suggests 10 as a 
typical flgvcce for a 40 module problem. ¥e may also 
put in the condition of predetermined niimber of iterations 
as stopping condition if such stability is not reached. 



39 


The program monitors line length during each 
iteration and uses a. simple averaging scheme of line 
lengtlTS from past iterations ( after a reorder) to compare 
the change in the present iteration. This removes the 
restriction that the process roust produce a non-increasing 
total line length and covers isolated cases of modxile 
oscillations. 

4.5 EXP AM) PROCESS 

This process has not been implemented in this 
project for two reasons^ that apart from complexity 
to the program <& time factor, the experiments with few 
circuits could not indicate whether an infinite loop 
situation will be avoided for all circuit problems. Secondly 
where an assignment to predetermined slots are to be made 
after this, the few circuits tried showed l,i 'tie difference 
to the cost matrix whether the assignment is done 
directly after converge or after expand just allowing 
■the modules to remove the overlap in converge. 

A brief description of this process follows. 

In this process the modules should assume their 
appropriate physical dimensions and their motion is 
restricted by boundaries of the board. The principal 
function is to increase the sizes of the modules and 
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reduce the physical overlap between all modules without 
destroying the orientati'in prepared by the CONyERG-E 
process. The paper suggests implementation of this 
process in two phases; 'bounded' and 'unbounded'. In 
the unbounded phase the modules are still not restricted 
by the boundaries of ci:rcuit board and physical overlap 
reduction is the main aim. In the next phase an additional 
vector force due to bounaary of board is included, 

4.5.1 REDUCING PHrSIO.AJi OVSEJilP 


For the modules to overlap the distance separating 
them in both X and Y sense must be less than half 
the sum of the dimensions in either the X a,nd Y sense. To 
separate this we depend on iiie vector component of 
repuD.sive . forces which a,re based on overlap in both X and Y 
sense. Eigijire ,4.7 below shows the physical overlap and 
the repulsion forces. 



I 


to ur>iiTS 



Eig. 4.7 Physical overlap and Repulsion vectors 
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If a pythogorean vector is used the movement 
would he effective. But uae of only the smeller of X or 
Y components (in this case X) of repulsion would provide 
two advantages; Firstly overlap could be reduced more 
rapidly. Secondly the resu.lting conditions would be 
more accurate approximation of the original layout (after 
converge) and it is important to preserve the relative 
positioning derived by COITl/lRG-B. Fig. 4.8 below 
il.lustrates this. Mhen they overlap equally resultant 
vector sum can be used. 


A 



Fig. 4.8 Use of Resultant vector and smaller of 
i and Y vector 

4.5.2 OTHER CC ITS I DERATIONS 

Once these 'physical' repulsion vectors'are computed 
for all modules, the 'electj^ical' attraction foi’ces (due 
to net interconnection) or vectors as- in converge are 
added to each module's total force vector. Certain constraints 
can be added to insure that overlap is reduced. It is 
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important to reorder nets between iterations but here 
we may use a simple constant (say 6) as reordering 
frequency because line length, changes very little 
between iterations. 

4.5.5 STOPPIUG CONDITION 

Since the reduction in overlap is the main aim, 
that, and not the line length change, f onus the criteria 
for stopping. Further, line length does not change much 
when certain value of minimum level of overlap is reached^ 
(for example when overlap is less than 20 percent for the 
modu3.es) the process then can be terminated. 

.4.6 FINAL ASSIGMENT 

The problem of ass igiii-jig the modules to sites that 
were originally defined in the input data is different type 
of distance minimizing prob3.em. At the end of converge/ 
expand each of the 'n' modules has a cost (s^j) assignment 
(simply the absolute value of rectilinear distance to the 
site) to each of the 'n' defined sites. Bach module has to 
be assigned to exact3.y one site and only one module is 
assigned to each site. If sites are extras dummy modules 
can be assumed with no cod: of assignment. A cost matrioc 
whose elements are a . . (l 1 it 3 £ is pi’esented to Ihe 
assignment algorithm for deciding appropriate assignment 
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( without, trying all the n.i cisihle assignments) such 
that is minimum where p is a particular 

permutation of first n integers that corresponds to the 
assignment of sites to modules. 


4.6.1 THE NlTURl: OP SOLUTION 

It should be noted 'Jhat instead of minimising 
interconnecting wirelength, the final assignment minimises the 
sum of the distances each module moves when assigned 
to a final site . In other words, the original criteria 
of electrical interconnecting lengths is ignored by this 
algorithm and hence its success depends on the approximation 
presented to it. That means the positioning at the end 
of converge /exp and should be a good approximation of the 
final answer, causing only a small, adjustment of movement, 

4.6.2 HUNG.^,RIM METHOD 


Basically the method successively modifies rows 
and columns, of cost matrix until there is at l.east one zero 
component in each row and each column such that complete 
assignment can be made against n independent zeros 
in the altered matrix. This will be an optimal assignment 
in the corresponding positions in the original cost 


matrix since on the basis of Nonig’s theorem a constant 
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can be added to and subtracted from argr row or any 
column withcu-^ changing the set of optimal assignments 
(The matrix elements must be non-negative intergers); 
and the theorem that if the minimum niunber of 'lines’ 
required to 'cover' each aero element at least once is 
'n' then there exist 'n' independent aeros in the 
matrix and vice versa. 

4.6.5 STEPS OF TEE HUIGEPvIAIT AttGORITHM 

¥e are given a cost matrix (elements a. . are 

i 3 

non-negative intergers re|)resenting the absolute value 
of rectilinear distance of a module i to site j). Through 
the following steps we can find the 'n' assignments 
that will minim iae the co-~t. 

Step 1 ; If the mini.m-um element in row i is not zero, 
then subtra-ct this minimum element from each element in 
row i (i = 1,2,...,m), 

Step 2 ; If the miiiim'um element in column j is not aero, 
then subtract this minimum element from each element in 
column 3 ( 3 = 1 ,2, . . . ,m) . 

Step 5 : Examine vovs successively, beginning with 
row 1 , for a row with exactly one unmarked zero. If at 
least one exists, mark this zero with the symbol. (<^X) to 
denote an assignment. Cress out (X) the other zeros in the 
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same column so additional assignments will not made to 
that coliimn (activity).. lepeat the process until each 
row has no unmarked zeros or at least two unmakred zeros. 

Step 4 : Examine the cblumns successively for single, 
unmarked zeros and mark them with the symbol (iX) to denote 
an assignment. Gross out (1) the other zeros in the 
same row so the corresponding resource will be assigned 
to other activities. Repeat the process until each 
col-umn has no vmmakred zeros or has at least two unmarked 
zeros. 

Step 5 ; Repeat steps 3 and 4 successively (if 
necessary) until one of three things occurs: 

(a) Every row has an assignment (A). 

(b) There are at least two unmarked zeros in each 
row and each colxrmn. 

(c) There are no zeros left unmarked and a complete 
assignment has not been made. 

Step 6 : If (a) occurs, the assignment is complete 
and it is an optimal assigimaent. If (b) occurs, arbitrarily 
make an assignment (A) to one of the zeros and cross out 
(I) all of the zeros in the same row and column, and then 
go to step 3. If (c) occui-s, go to step 7, 
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S tep 7 ; Check (v^) all rows for which assignments 
have not been made. 

Step 8 ; Check columns not already checked X'rhich 
have a zero in checked rows. 

Step 9 : Check rows not already checked which have 
assignments in checked columns. 

Step 10: .Repeat steps 8 and 9 until the chain 
of checkings ends. 

Step 1 1 ; Draw lines through all unchecked rows 
and through all checked columns. This will necessarily 
give the minimum number of .lines needed to cover each 
zero at least one time. 

Step t2 : Examine "he elements ti^t do not have at 
least one line through them. Select the smallest of these 
and subtract it from every element in each row t.hat 
contains at least one uncovered element. Add the same 
element to every element in each column that has a vertical 
line through it. Return to step 3. 



CHAPTER 5 


PLACEMENT IMPLEMENTATION 

The force-vector heuristic technique discussed in 
Chapter 4 has been implemented in Pascal (Appendix II). 

The following is documented for easy understanding of 
the user. 

This chapter will discuss data input, program 
overview and data output in the next three sections. Before 
discussing data input, it should be noted that the board 
dimensions are program constants and hence should be decided 
first and changed accordingly for each circuit. The 
difficulty in reading them from input file is because of 
'window ‘and 'view port' dimensions which will also be 
program constants depending upon the Board dimension chosen. 
Therefore before executing the placement program for each 
circuit the board dimensions in 1 and Y directions and the 
window and view port dimensions are changed through 
the use of a Mic Pile, (Use of Mic file discussed in sec.to**f' 

Depending upon the circuit and routing requirement 
the Z-span and Y-span of board are chosen as some integer 
units. Each integer unit may be typically or 1 or 2 mm' 
There is no restriction on the size (excepting window and 
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view porting constraints) as far as the placement problem 
is concerned* But if router program is to be used then 
the dimension is restricted to within 140x100 integer 
units. 

Apart from these program constants of X~span, 

Y-span and the corresponding window and view port dimensions 
the following data are also to be supplied for each circuit 
for placement program. 

5.1 DATA INPUT 

This can be divided in three parts in the same 
sequence as follows. All measurements of length or size 
are interms of the integer units decided for the board 
dimensions. Appendix I contains a listing of input for 
circuit 1 , 

Part 1 ~ Module data ; The sequence of data under 
this part are: 

1 , Number of modules 

2. Individual module data (modules are numebered 
in the order in which these sequence of data are 
read for each): 

(a) !H5rpe of module (integer coding): 

1 : Plxed 

2;Ih*ee 

3:X PiAed 



4;Y Fixed 
5 jX Bar 
6:y Bar 

Size in X-direction and then size in Y-direction 
from the coordinate position (centre) of the 
module. Normal coordinate conventions used. 

Coordinate position in X and then in Y. Note 
that normal coordinate convention, i.e,, origin 
at left bottom comer of board, is used. 

Number of pins in the module (dual in-line). 
Orientation of chip: 

0 : If not applicable 

1 : First pin on left top 

2 : First pin on right bottom 

The data in (d) above and this are not 
applicable for bar modules (and hence =0), but 
applicable to subcomponents of the bar module. 

The subcomponents are identical. In 
circuit 1 example the gro-und and power lines are 
modelled as bar modules and have no subcomponents 
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The data in (e) is not applicable for 
connector segment modules (and hence =0), But 
data in (d) is applicable depending upon how 
many board terminal pins are modelled as one 
(fixed) connector segment module, 

(f) Humber of submodulesS 

Applicable only for those bar modules 
that contain subcomponents. (Hence =0) for all 
other modules), 

When this data is applicable, this must 
be immediately followed by data in (b ) through 
(e) for each of the subcomponents. 

Part 2 - Het Data ; Hie seq.uence is as follows. 

1 , Humber of Hets 

2. Individual net data: (Hets are numbered in the 
order in which this data is read). 

(a) Size of net (i.e. number of pins in net) 

(b) Wiring rule for the net (integer code); 

t ; Star rule 

2 : Serial 

3 : Span (minimum spanning tree) 
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(c) For each pin in this net: 

i) Module number 

ii) Pin number 

iii) Sub-module number 

Sub-module number is applicable only if 
the connection is to a subcomponent of a bar 
module whose number is given in (i). For all 
other cases it will be zero. 

Note that the placement program does not 
move submodules. In this program submodule 
pin numbers are assigned before-hand to a net. 
They can be easily modified to be specified at 
the end of the placement. In any case the 
actxial location of the pin is required for 
routing. 

As far as the placement program is 
concerned, for calculation of vector forces in 
this net only the perpendicular distance to 
the concerned bar module will be taFen into 
account. 

It is important to ensure that, first, 
the nets read are nets of power line and 
ground line . 
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Part 3 -■ Details for Terminal pins of Board ; The 
sequence under this is as follows; 

1. Total number of terminal pins of board, 

2, Length of a terminal pin. (The scale for this 
is the same as module size or position 
coordinates), 

5, Orientation of Board terminals (integer code); 

1 : Terminals on top of board 

2 ; Terminals to right of board. 

It should be noted that if we are to use the router 
program with orientation 2 then it requires modifi- 
cation to one procedure (Power and Ground lines) in roster 
program, which otherwise assumes orientation 1, 

For placement program alone both can be used. This 
helps to adjust the window parameters so that the 
longer side is along X-axis. 

4. X coordinate and T coordinate positions (bottom 
point) of all terminal pins in the order which 
they are numbered on board, 

5.2 BATA OUTPUT: (Eefer Appendix III and Appendix Till) 

One part of this is tailored to use the router program. 
Ihe board with the given dimension is thought of as a grid 
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of cells. Ihe number of cells in the grid are equal to 
(X span*T span),. As mentioned in the beginning of this 
chapter it is restricted to (140x100) if router program is 
to be used, 

.Another part of the output is to produce the final 
placement configuration on the screen or plotter. This is 
done by procedure ’Draw Final Placement’, 

Yet another part pertains to output statements which 
have been suppressed with comment delimiters. These are 
obvious on inspection of the program and can be made to 
output on a separate file. These are meant for talcing a 
print out of various forces on modules, their resultant etc, 
and the cost matrix formed for assignment problem and 
the alterations in the matrix by the assignment algorithm, 
only that part of output pertaining to router program is 
discussed further. This in turn can be thought of as the 
following three parts in the same sequence. 

Part 1 - Board Dimensions ; The sequence is; 

1. X-^span of the board 

2, Y~span of the board 

These pertain to the board dimensions discussed in 
the very beginning of this chapter and are constants of 
placement program for each circuit. They are merely 
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presented at the output again for use "by router program, 
routing for the same circuit. 

Part 2 - Cell Uumhers of all pins : The sequence is: 

1, Sequence of cell numters corresponding to all 
pin locations of all modules (including 
subcomponent pin locations of bar modxiLes). 

2, Integer 99998 to indicate that part 1 and part 2 
(of input to router program) are terminated. 

The formula for knowing the cell number of a pin 
location is ( (Y coordinate position of pin-l)*! span) + 

(X coordinate position of pin). The cell niunber will be 
an integer since centre of a module is an integer position 
and its size is chosen as an integer depending upon the 
ntuaber of pins. Therefore each pin location is an integer 
position in X and Y. It is important to note that the 
choice of coordinate position in placement program makes the 
left bottom most cell as cell Uo, 1 . But the router program 
has left top most cell as first cell increasing linearly. 
Therefore a conversion is required in coordinate position 
of pin before giving the correct cell number to router. 

This is automatically done in placement program. 

The order in which the sequence of cell numbers are 
read is unimportant since they are required by router only 
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for the sa]£e of marking these cells as occupied by pins 
and to physically define a module. However for convenience 
an order is adopted in this output inplacement. 

Modules are examined one by one in the order in 
which they were input and all pin locations for this 
are output. When a bar module is encountered, if it has 
submodules then all the pin cells of all its subcomponents 
are outputj if it has no subcomponent (for example, a power 
or ground line) then no cell number is output). 

As soon as the first connector segment module is 
encountered the cell numbers of all board terminals are 
output. 

Part 3 ~ Cell numbers of Pins of each net : The 
sequence is; 

1 . Sequence of cell numbers to be connected to all 
pins of power line followed by integer 0 to indicate 
the termination of cells in this net (net l). 

2. Similar sequence for let 2 (ground line) 
followed by integer 0, 

5. Similar sequence for let 5, let 4 and soon till the 
last net, Bach sequence is followed by integer 0. 

, Last net sequence is followed by integer 0 and then 
integer 99999 to indicate end of data under this 
part. 


4 
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The order of the cell numbers in each sequence for a 
net depends upon the order in which they are stored in the 
linked list for that net in the placement program. Both this 
order and the order in which the net (i.e. the sequence 
of cells) are considered for routing may he changed hy 
router program. 

Changing the order within a sequence of cell numbers 
corresponds to implementing a wiring rule. Placement 
program can he modified to do that after final placement 
and output cell numbers. But that would alter the router 
program. The router implements some wiring r\ile by itself 
and a few more can be added to it, 

5.5 PROGRM OVERVIEW 

5.3.1 STEPS OP THE PROCESS 

Step 1 : Initialise the device, window and view port. 

The window and view port dimensions are fixed by' the ’init 
procedxare'. 

Step 2 : This is done by procedure 'Readin', in the 
following sequence : 

(a) Read module data, from input file (Ref. Sec, 5.1). 

The details of each module is stored in the variable 
arra.y 'module* including Bar Modules as 'component' 
records. If the bar module has any subcomponent 
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the details of all its subcomponent are stored in 
the variable array ’Barmap* as *Bar component* 
records. 

(b) Fix the average size of the modules.. Siis is 
required for calculating repulsion force:. 

(c) The candidate modules and candidate sitep are 
fixed for the assignment problem. The candidate 
modules are of free type or X-fixed or Y-fixed 
type. These module numbers are stored in the 
variable array ’Matmods* and their X and Y 
positions are stored in the variable array 'Mat- 
sites',- i.e,, the candidate sites. The Recti- 
linear distance of each candidate module to each 
of the candidate sites as at the end of placement 
problem decides the assignment matrix. 

(d) Read Net data from input file. (Ref. Sec. 5.1). 

For each of the net read in, form a linked list. 
Other details of the net are stored in the 
variable array 'net' as a 'net-cell*. Each net- cell 
has a pointer to the appropriate' linked-list. 

(e) Note down the number of pins in the net with 

4 

maximxmi pins against variable 'maxpin net'. This 
is used for progressive net addition. 
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(f) Read the terminal pins data (Ref. Sec. 5*1 ). Store 
the X and Y positions of each of tie terminals 
pins of the board in the variable array 
’ Termpinlocs* , 

atep 5 : Implement the procedure ’converge’ discussed 
in detail in Chapter 4. This is done by procedure ’converge* 
which calculates the vector forces on each module and moves 
them as discussed in Chapter 4, by mainly using its two 
local procedures ’Net Reorder’ and ’Move’. The procedure 
’Net Reorder’ implements the specified wiring rule for 
each net by deciding which module (pin) is to be connected 
to which module (pin). It then updates an important array 
variable ’Temp' that stores ih its records all the connection 
details. The procedure ’Move’ then calculates all the 
vectors on each module using ' Temp * variable and also the 
resultant and mobility factor in Z and I directions separately. 
Next it adds the change in position to the appropriate field 
in each component record in the array ’Module’. ’Converge’ 
procedure is illustrated in Rig. 5.1, 

Step 4 ; 3he assignment process discussed in detail in 
Chapter 4 is implemented by the procedure ’assignment’ 
mainly through two procedures local to, it: The procedure 
’Fix-matrix’ forms the cost matrix for the assignment 
problem (Ref. Step 2(c)). It should be noted that in the 
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Fig. 5.1: Plow' chart for Step 3 (converge) 
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case of X-fixed (and Y-fixed) modules the cost of assigning 
it to any site not haring the same X-coordinate (I-coordinate) 
is made arbitrarily large, !Ehe procedure 'Hungarian' then 
takes orer and implements the steps 1 through step 12 of the 
Hungarian assignment algorithm discussed in Section 2.3 
by means of a number of procedures local to it- 

Step 5 ; Draw the final placement configuration on the 
screen/plotter. This is done by procedure 'Draw placement' 
which invokes it two local procedures 'Mark pins' and 
•Mark termls' and also to procedure 'Draw Reef. All tiK 
modules (except subcomponents of Bar modules and connecting 
segments) are drawn in their final positions on the 
board and numbered. The pins of modules and the terminal 
pins of the board are drawn and their first pin is 
numbered. 

Step 6 ; Output the program constants X-span and 
Y-span of the board for use by router program. 

Step 7 : Examine mod Tiles one by one. Output cell 
numbers of all pins (Ref, Sec. 5,2) of all modules considering 
) the first cell to be on left top comer. This is 
used by router program. No pin-cell number output for 
bar modules with no subcomponent. \i/hen the first connector 
segment is encountered the cell numbers of all terminal 
pins is output by accessing the variable 'Term pin Iocs'. 
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This step is implemented by procedure 'Write all pin 
cells ' . 

gtep 8 ; !I3iis is the last step in placement program, 
implemented by procedure 'write net pin cells’. Examine 
nets one by one and output cell numbers of all pins in 
each net (Ref, Sec. 5.2). 

5.3.2 SOME IiyiPORTANT TlRIiBLES AND GORSTAITS 

Avlength ; Global. Gives the average value of total 
line length fc-e: v; moves (excluding the present move) 

since the last net reorder. Ref. Section 5.3.1, Step 5 
flow-chart, 

Avsize : Global. Gives average size of modules. Used 
for calculating reptxlsion force. Ref. Section 5.3.1 Step 2(e) 

■ AVX(MX and MOBPACTX ); Real variables local to 
procedure 'Move*. Once the resultant force on each module 
is calculated, the average force on modules, maximum force 
on modules and the mobility factor on modules are calculated 
against these variables (in X direction). Similarly for 
AW, Ml and Mob Pacty. 

Bar Map (Bar comp and Bar ind) ; Global. Illustrated 
below. The upper bound of 20 for the array can be varied 
as required. 
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Pig. 5.2 

_ bxT : Local to procedure 'force*. A Boolean array 
to indicate whether a module is connected or not to each of 
the other modilLes. Updated every time a particular module 
is considered. 

PiffFY (ORIEHT and PIULEN): Global, but used only in 
procedures 'Markpins* and 'Write-cells'. Pertains to 
distance between pins in a module and should be assigned a 
proper value in these procedures by the user. 'ORIENT’ 
refers to orientation of chip (Ref. 5.1 Part 1(e)). ’PINLEN’ 
is a constant local to procedure 'Mark pins' to decide the 
length of pins. User specified. 

FX( and PI) t Array variable local to procedure 'force 
and used for storing all the forces on a module updated for 
each module. This array is passed on as parameter to 
function resultant which calculates the resultant force in 
X-direction (and Y-direction) on the module by sorting them 
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and deleting pairs of opposite sign as discussed in 
Chapter 4. I!he same identifiers are used as real variables 
local to procedure move, in which case, it is merely used 
as a temporary storage to aid in calculation of maximum 
force in X direction (and Y direction) on modules, 

Preq , ; G-lobal variable whose value is decided by 
procedure ’compare’. This integer value is the freq.uency 
of iterations of converge between the net reordering. 

From , UP TO ( and Map ) ; Array variables local to 
procedure ’Net reorder'. ’From’ and 'Upto' give indices 
to ’Map’ which stores module numbers which are to be 
connected in a nit updated for each net considered, 

INC : Global, Boolean variable to decide whether a 
particular net is to be included for consideration during 
progressive net addition of the ’converge' process. 

Matrix. Mat Mods and Mat sites : Matmods and 
matsites are global arrays discussed in Sec, 5,3.1, step (c) 
^ 'Matrix' is the cost matrix formed from these for the 
assignment problem. A row of the matrix pertains to 
distance of a module to all the sites and other information 
(whether the element is unmarked crossed or assigned 
during the process of assignment i.e, 'Daagtyp'). Also 
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TEBMl TYP ; Specifies orientation of board terminals. 
Refer Section 5.1» Part 3. 

TEMP ; In important arraj variable that stores all 
interconnection data between modules. ’Tempmax* can be 
changed by user. Each module has a pointer to one cell 
of ‘Temp' as illustrated below. 



TEBMLOCS ; Array storing positions of all board 
terminals. 

5.3.5 BRIEF ROTES ON PROCEDURES 

Most of the procedures have been outlined in 
Section 5.3.1. A further, note is added here. 



static Uesting of Procedures and Punctions 


MAIN 

DHAWBECI 

CONVERGE 

REAEIN 

REAEMODULE 
PIXAVSI ZB 

FIX MATMOES AND SITES 
READNETS 
FIXMMPIN 
REAETERMLS 
DISTANCE 
LINELENGTH 
NETREOEDER 

STAREEORDER 
SERI AD REORDER 
SRMREORDER 
INSERT 
MOVE 

FORCE 

RESULTANT 

SIGN 

XVECTOR 

YVECTOR 

XREPULSE 

YREPULSE 

ROUND 

COMPARE 

ASSIGNMENT 

FIXMATRIX 

HUNGARIAN.- 

FIXSTATUS 

MA.TPRINT 

STA.TRPRINT 

STATCPRINT 

ZEROPRINT 

ROWMIN 

ASSROW 

ASSCOI 

CONDITIONCHECIC 
REDUCE MATRIX 
TICKLINES 
COVERLINES 
DO FINAL MAT 
ASSIGN 
WRITE CELL 
WRITE BAR CELL 
DRAW PLACEMENT 
MARK PINS 
MARK TERMLS 
WRITE ALL PIN CELLS 
WRITE NET PIN CELLS 


Fi«* Ne45ting of Procedures. 
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Procedure Fet Reorder s Refer Sec. 5.3.1, Step 3. This 
implements thi^ee wiring rules. In ’star' wiring rxlLe ihe 
module in the central position of the net (or a Bar module 
without suh components ) is connected to all other modules in 
the net. 'Span' Reorder implements a minimum spanning tree. 
In Serial Reorder every vertex has an incidence of utmost 
2. In all, if a net ha^s 'N* pins then there are (R-l) edges 
to be specified using 'Prom' and 'Upto' array. At the end 
of net reorder the procedure 'Insert' updates the array 
variable 'Temp'. 

Procedure Move ; Examines each module. Prom 'Temp* 

array information calculates attra,ction force to each 

connected module and repulsion force to each unconnected 

module specified by ’Bit' array. The functions X vector, 

Y vector^ X repulse and Y repulse calculate these forces 

with due regard to module types. The local procedure 

•force' uses these functions and the function 'resultant' 

to update the Force field on each module. Once this is 
done 'Move' calculates the mobility factor and 
updates the position field of individual modules. 

Pun ct ion Distance ; Given two modules as parameter, 
this calculates the rectilinear distance between them with 
due regard to module types. 
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Procedure line length ; Uses ribove f-uinction rnd 
connection details in 'Temp* array. Returns the total line 
length for the current lay-out a,nd net reorder imple- 
mentation. 

Procedure compare ; Compares line length after every 
move with the average line length in previous moves since last 
net reorder. Depending upon the absolute value of the percent 
change in this compf^rison it decides to stop the iteration 
or the frequency of iterations for the next reorder. 

Some Local procedures of Assig:nment ; The procedure 
'Fix Matrix'^ forms the cost matrix for the assignment prohlem, 
'Fix status* initialises the matrix. The 'print' procedures 
are added for diagonostic purposes and have been 
suppressed. 

Procedure Hungarian ; This implements the assignment 
algorithm of Section 4.6.5. The procedures 'Row Min' 

'colmin' correspond to steps 1 and 2 of algorithm. The 
'while' loop containing procedures 'Ass Row', 'Ass Col' and 
'Condition Check' corre^ond to steps 5,4 and 5 and 6 of 
algorithm, \lhxle within the loop it is analogous to out 
come in step if condition check causes exist from 

while loop due to success it means the out come in step 
5(c). If the cut come oorreeponding to step 5(c) is caused 
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then the procedxire 'Reduce Matrix' takes over, which 
corresponds to step 7. The repeat loop in this case 
ensures that after step 12, we return to step 3 of 
algorithm. The inner procedures represent steps 8 to 12 
of the algorithm. The first repeat loop within 'tick lines' 
implements steps 8-10 in the form of two inner 'for' 
loops and an inner 'repeat' loop respectively. Procedures 
'coverlines' and 'Do final mat' implement steps 11 and 12 
of algorithm. 

Procedure Assign; In case of success this procedure 
updates the position field of individual modules. 

Procedure Draw Rect t Draws all the modules except 
subcomponents of Bar modules called by converge process 
to give a visual pictxare of converge to the user. Procedure 
draw placement also makes a call to this. 

Procedures to Output Cell numbers ; It should be 
noted that connector segment moduLes (i. e. , those containing 
board terminals) are identified by their 'orient' field * 0 
but pins <>0 where as bar modules are identified by 
orient = 0 and pins =0. In case a bar module has 
subcomponents its barsubs <> 0 and this field is zero fDr 
all other moduJ.es. The procedure 'write cesll ' outputs 
a particular cell number when the module number and pin 
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number (may be a board terminal nmber) are specified. If 
pin number is specified as zero in iiie case of a connector 
segment then it outputs cell numbers of all the board 
terminals, procedure ’write bar-cell' outputs the 

cell number of a particular pin of a subcomponent when the 
bar module, its subcomponent number and pin of subcomponent 
are specified. If the last two parameters are zero then it 
outputs cell numbers of pins of all the subcomponents of that 
bar module, These two procedures are invoked by 'widte 
all pin cells' and ’write net pin cells’ to produce output 
data as discussed in Section 5,2. 

5.4 USE OF 'MIC FILE FOR FRO&RAM EXECUTION 

As seen in Section 5. It for each circuit the program 
constants of window and view port dimensions are to be 
changed depending upon board dimensions. The board dimen- 
sions are also kept as program constants for this reason. 

Further when the placement program router program 
and the plotting routine are to be done in sequence or 
independently and that for many circuits, extensive renaialng 
of files will be required to avoid input output files 
confusion. For this reason as soon as the program is 
executed the input file is stored back as PLAl.IN (i.e. 
for circuit U circuit 2 it is PLA2.IN etc) and the 

placement output as PLA2.'OUT etc. 
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These require a series of instructions for 
execution and are convenient to do through the use of a 
'MIC’ file. For each circuit a Mic file instruction is 
necessary. For example, the 'P lace 1, MIG' file instructions 
require a shnple command 'DO PD1CE1' to carry out necessary 
instruction for placement problem of circuit 1. 'Place2' 
for circuit 2 etc.), 

'Ehis can be immediately followed by MIC file 
instructions * Route 1 ' and 'Plotl ’ and then 'Ramel* 
instruction which renames the input/output files as they 
should exist. If the user wishes to terminate his work 
with placement or in any other combination 'Ramel ' must 
follow 'Placel*. 

The Mic file instructions for circuit 1 are 
enclosed as Appendix VII. 



CHAPTER 6 


ROUTER IMPLEMEIMIOR 

The router presented here (Appendix IV) is an 
existing impleiaentation [10], This chapter merely attempts 
to document it more clearly for ea.se of user and to piojvide 
an environment integrated with the placement program. 

6.1 DATA IIPUT 

Ttiis can be thought of as a sequence of four parts. 

The first three parts and the sequence within each part 
is identical to that of placement prcgran output. Since 
the input data can be given independent of placement, the 
sequence is repeated again and points pertaining to router 
program alone are briefly discussed. 

Part 1 Board dimensions : 

1 , X span of board 

2. T span of board 

These are read against router program variables 'X Cooird’ and 
*Y Coord*. Program considers the board as a grid miii 
(X span * Y span) cells limited to 140x100. This can be 
changed depending on memory available. A cell is a 
sqtiare or 1 or 2 mm square depending on circuit design 
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typically. The first cell is assumed at top-left corner 
increasing towards ri^t and then on second line and so on. 
The last cell is at bottom-right corner. 

Part 2 ; The sequence again is; 

1. Cell numbers of all pins of all modules 

2, Integer 99998. 

The order in which the cell numbers are read is unimportant 
as they only physically define a module. The router 
program recognises only cell numbers without any regard to 
type of module. 

Part 3 - Cell numbers .of pins of each net ; !Ehe 
sequence is: 

1. Cells of all pins connected to groimd line 
followed by integer zero. 

2. Cells of all pins connected to power line followed 
by integer zero. 

5. Cells of pins in each of all the other nets. Each 
net is a sequence of cell numbers and followed 
by integer zero. 

4. Integer 99999. 

The ordar in which each sequence under (3) above is 
considered may be changed by the router program by means 
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of either of two procedures ’Area ordering* and ‘Mid ordering*. 
The order in which pairs of cells are considered within 
each sequence may also he changed by router program by- 
means of function 'nearest*, which infact is an implementa- 
tion of minimum spanning tree wiring rule. 

The ETets 1 and 2 are implemented as a *star* 
wiring rule in which each of the pins is connected to the 
ground line or power line as the case may be. 

Part 4 - Character codes as data ; The program 
expects the following sequence of single characters as 
code for choosing methods of routings 

Character *M* or any other character s If *M* is 
read the nets are ordered for" routing depending upon 
which net (any pin of it, i.e. any cell within 
the sequence of cell numbers) is farthest from 
centre of board by invoking procedure *Mid oi^er*. 
If M- is not read then the nets are ordered for 
routing depending upon which covers most area of 
board by invoking procedure ’area ordering*. 

2. Character *B * or any other character ; If ’B* is 
read then the ordered nets are considered in 
descending order (of distance or area covered) for 
routing. If ’B ' is not read then they are 
considered in ascending order. 
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5. Character *C* or any other character ; If ’C* is 
read then 'chaining method’ is used for finding 
a path between given pairs of pins. If ’C’ is 
not read then 'Steiner’s tree’ (i.e* procedure 
’Straner's tree' ) is used for the same purpose. 

4. Character *0 * or anv other character ; If ’0* is 
read then the program does not route for the 
other side. If 'O' is not read routing is done 
on both sides. 

6.2 MIA OUTPUT 

The 3X!uter program operates on the input data to 
produce the following data strictly in the same order. The 
actual drawing on the screen or plotter is not done by 
the router program but by a small plotting routine in fortran 
which is discussed in Section 6,5. This data again C 0 ,n be 
considered in four parts in sequence. An example is 
appended. (Appendices V and 1/1), 

fart 1 •» Data to Ground and Power Lines ; Por each 
line 5 points (cells) are output including the corner 
point. Though line drawn covers only one cell in width, 
it is actually two cells in width. The sequence iss 

t. Three cell numbers which are to be joined by two 


straight lines in the same sequence followed by 
integer 0 as delimiter. 
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2, iDiree cell numbers similar to above, but which 
make the power line. 

5* Integer 99996 to terminate part 1 data. 

Part 2 - Routing I>ata for side one ; Ibis mainly bas 
many sequences of cell numbers. Sae cell numbers in each, 
sequence are to be connected by straight lines in the same 
order by plotting routine to make the routing. Each cell 
number in one sequence does not necessarily denote a pin 
location, but at least one and at most two of these cells 
must correspond to pin locations, since these sequences are 
formed by either of two procedures ( ’ chaining method ' or 
‘straner’s tree*) which receive a cell pair as their 
parameter. Hence sequence iss 

1, Sequences of cell numbers, each sequence 
followed by integer 0 as delimi.ter. 

2. Integer 99997 to terminate data in this part. 

Part 5 -• Routing Data for side two : This is . similar 
to part 2, But produced only if desired by user through 
the character code at input as discussed in Section 6.1. 
This part, whether output or not, is terminated by 
integer 99999. 
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Part 4 - Data of -onmade connections : This part is 
generated on a separate file called ’out' since it does 
not concern the plotting routine. This lists out in 
sequence for each unmade connections: 

1 . Net number 

2. Pin number in the net 

These connections have to be manually made, in 
example of router output is included as Appendix V. 

6,3 PLOTTING ROUTINE: (See Appendix VI and IX) 

Though the router program produces a sequence of 
cell n'umbers to be actually connected, it does not do the 
actual plotting or drawing. This is done by a small 
Fortran Routine using just 'line' cells of GPGS. 

This routine uses two data files which are nothing 
but the output file of placement program i.e., same as input 
file of router program and the output file of router program. 
These two are renamed FOR 20, DAT and FOR 21.DAT respectively. 

FOR 20, DAT is only partially used, i.e. to mark 
short dashes at cell positions corresponding to all pins 
on the circuit board, i.e, the cell numbers after integer 
99998 in the file are not used. 

FOE 21, DAT is used as follows: 
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(a) For drawing power and ground lines, i.e. cell 
numbers upto integer 99996. 

(b) For drawing all possible wire connections on 
first side, i.e. data between integers 99996 
and 99997. 

(c) For drawing interconnection on the second side, 
i.e, data between integer 99997 and 99999. 

Before drawing connections on both sides the program 
awaits for user response (type any integer) and also for 
clearing device/plotter. 

6.4- PROGRAM 07ER VIEW 

6.4.1 STEPS OF THE PROCESS 

Step 1 ; Read X span and Y span of the board. 

Step 2 ; Read cell numbers of all pins of all modules 
from input file (Ref, Sec, 6,1, Part 2), Update the 
information for these cells as follows: side one occupied 
and the cell not beii^ used ty any particular path. Say, 
the fields for each cell used for this pu 2 :p)Ose are; 

Cell Info [Cell number]. side = 1 
Cell Info [Cell number]. Path code =1 

One more field i.e., cell info [cell number], Link cell is 
not of concern at this step. 
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Step 3 : Read cell numbers of pins of all nets {Ref, 

Sec. 6.1, part 5). Form a ’net matrix*. Each row represents 
one net and each element of this matrix ha.s two fields: Cell 
number corresponding to the pin and Boolean info whether 
the cell has been tried for routing by 1he program (i.e. 

0 or 1 for 'irntried* and ’tried’ respectively). Store the 
cell numbers of pins in the net in the appropriate element 
and update the other field to zero (i.e. Efot yet tried for 
routing. Ihe elements of this net matrix are accessed as 
net matrix [Uet ho., Pin No.] where net number is the order 
in which the nets were readin and pin number is the order 
in which the pin (i.e., its cell number) was readin. The 
end of pins in a row is indicated by a zero stored in the 
cell number field of the element. 

Step 4 ; Note down the total number of nets. 

Step 5 : For all border cells update cell info as in 
step 2. This is done by procedure ‘initialization’ . 

Step 6 : Demarcate power and ground lines. All 
second and third cells from left boundary and second and 
third cells below top boundary upto a distance (in this 
case 16 cells from mid point) are encoded as in step 2 
but here the PATH code = 2, indicating that they are to be 
used for connections of power net (i.e. actual net number one). 
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Similarly ground line uses second and third cells at right 
and below top boimdary with PATH code = 3. This step is 
implemented by procedure 'POLIHE ' . 

Step 7 : Output three cell numbers (i4e,,end cells 
and corner cell) of power line followed by integer 0 and 
then three cell numbers of ground line followed by 
integer 99996. 

Step Perform the following steps for side one 
(Program Variable ’Side’ = True). 

Step 9 ; Consider ‘PATH’ Two; i.e,, connections to 
power line. Examine each element of first row of net matrix. 
Treat each cell number as ’start cell’ Eoid invoke procedure 
’Straners tree' with 'Start cell' and 'end cell’ as para- 
meters every time. Here the 'end cell’ can be taken as 
2 for parameter's sake. Ho need to change end cell’s path 
code since this is already fixed in procedure POIIHE. 
’Straner's tree' implements Lee's Algorithm and finds a 
path between the start cell specified and any of those cells 
with the same PATH code as Path under consideration. When 
a Path is found it outputs a series of cell numbers to be 
interconnected by straight line for this path (Refer 
Sec. 6. 1 ; part 2) , 
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Step 10 » Consider ’Path’ Three. This is siTnilax to 
step 9 but connections to groxmd line are made examining 
the second row of matrix (i.e. actual net number two), with 
3 as the path code for cells and as end cell parameter. 03ie 
output in this step is similar to step 9. Steps 9 and 10 
are in effect implementation of star wiring rule for the 
first two nets. 

Step 1 1 ; Read character code and order the remaining 
nets (i.e.. net 3 onwards). By invoking procedure ‘Mid- 
ordering’ if character ’M’ is read or procedure ’Area ordering' 
otherwise. In the case of ’Mid ordering’ the nets are 
ordered depending upon the maximum (Rectilinear) distance 
(for any pin of a net) covered by 1he net from the mid point 
of the board. In the case of 'Area ordering’ the nets are 
ordered on the basis of board area covered by each net (i.e. 
greatest X span covered * greates T span covered by the net). 
The net number and its associated area or distance is 
stored in the ordered sequence in a ’Met order array 

Step 12 ; Read two more character codes. If the 
first character read is *B ’ then the 'net order array* will 
be examined in descending order. Otherwise in ascending 
order after sorting. If the second character read is *C' 
then the procedure ’chaining method * will be used for finding 
path between cell pairs; otherwise Steiner’s tree method 
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will be used. ’Chaining method* is similar to ’Straners 
tree* except that the path has to end at the cell number 
specified as ’end cell’ instea-d of at any cell having a 
same ’path code’ as the path under consideration. Both 
req.uire ’start cell* and ’end cell’ parameters with the 
path nxmiber (all connections in one net have same path 
number) specified. While routing on first side, both store 
start cell, end cell and path nimiber information in a 
’second side’ array (along with the net number and pin 
number corresponding to start cell) if a path can not be 
found. 

Step 13 : Repeat the following for each net in the 
’net matrix* in the order decided by the ’net order array’ 
which yields a ’row number’ (net number) of the net matrix: 

(a) In the selected row, choose the cell number of the 
first element (i.e., first pin of this net) as 
start cell. Updated the ’tried’ field of this 
element to 1 , 

(b) Store the pin number of the tried pin in this net 
in a ’pin nimiber array’ sequentially. 

(c) Repeat the following until all pins of this net 
have been tried: 
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i) Implement a minimum spanning tree starting 
with this first pin. This is done hy 
function ’nearest* using a ‘local array’, 

ii) Each time a new vertex is added to 1±ie 

spanning tree take that cell number in the 
net as 'starting cell' and find the pin 
number (from 'pin nixmber array’ which stores pins 
already tried for routing ) nearest to it 

and take that cell as 'End cell', 

iii) Store the start cell’s pin number in the 
'pin number array' and update the 'tried' 
field in the corresponding element to 1 , 

iv) Invoke 'Straners tree' or 'chaining method* 
with these two parameters. 

The program assumes path number 5 for the first net routed 
in step 13 and increments this path number for every 
subsequent net routing tried. (Path number 4 also can be 
assumed in the beginning). 

Step 14 ; Perform the following steps for side two. 
(Program variable 'side' = false); output the integer 99997 
and initialise the 'second side' array index' to l . 
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Step 15 ; Read character code. If the character read 
is '0 ' then no routing on second side, else proceed 
with routing on second side. If no routing, then invoke 
procedure 'Efo lay' with 'second side array pointer' as 
parameter, incrementing it by three after every invocation 
since every third element of this array contains infomation 
on the start cell that was not connected (if any). If 
routing is to be done then invoke 'Strsner's tree* or 
'Chaining method' with 'Start cell' and 'End cell' 
infonnation contained in this array (if any) as parameters 
till such pairs are exhausted. At this time if 'Straner's 
tree' or 'Chaining method' is not able to find a path it 
will invoke proceduie 'No lay'. The procedure 'No lay* 
writes out (or 'out' file) the start cell, i.e., the pin 
numbers and net numbers which has to be manually connected. 
Pinally output the integer 99999. 

6.4.2 SOME niPORTANT VARIABLES AMR CORSTAFTS 

X Coord and Y Coord : These are the board dimensions. 
Ref. Sec, 6,1, part 1, Limited to 140x100. 

Array' 'A '; This is the 'net matrix* array discussed 
in step 3 of Sec. 6,4.1. The indices of the array can be 
changed by user to suit a circuit. 25ie 'select* and 
'data' fields of the elements of this matrix correspond to 
the 'tried' and 'cell number* fields in our discussion. 
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AC « An integer var parameter used in proceedures 
'Straner's tree* and ’Chaining method’ and indexes the 
element (in a list used by these procedures) that contains No, of 
the target cell (i.e,, ’end cell’ or a cell having the 
desired ’path code’). 

B: This corresponds to ’net order array' discussed in 

step 11 of Sec. 6.4.1. The fields of the elements of this 

array are self explanatory except the 'order' field, which 

corresponds to ’net number’ field in our discussion, 

one more 

There is also array B local to fmction 'nearest* 

mentioned as ’local array* in step 13 (c)(i) of Sec. 6,4.1. 

This is further explained under notes on procedures in the 
next subsection. The index bounds can be changed by user. 

Cell Status ; Corresponds to 'cell info' array 
mentioned in step 2 of Sec. 6,4.1. The fields 'select', 

'pin' and 'chain coord' respectively of each of its 
elements corre^ond to the fields 'side',^ 'path code* and 
'link cell’ mentioned in the same step. Index upper bound is 
limited to 14000. 

I Sig; This is an index to the ’net order array’ 

(i.e. array ’B’) that gives the net number to be 
considered for routing. 

jPOg ; This is an index to get the pin number (in the 
selected net) to be considered for routing as start cell. 
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LIS[PL t An array for use by 'Straner's tree* or 
chaining method to store all cell numbers available for 
laying a path between the start and end cells, (or the 
target cell). Index bounds can be changed by user. 

Sec Side : This index bound also can be changed by 
user. Array meant to store the start cell (&corresponding 
net ntimber and pin number), end cell and path number when 
the procedures fail to find a path on the first side 
between start and end cells. 

THBE s This is the ('tried* or connected) 'pin number 
array* pertaining to the net under consideration. Discussed 
in step 15(c) of Sec. 6.4.1, 


6.4.3 BRIEF NOTES ON PROCEDURES 

Static Nesting Sequence of Procedures : 

MAIN 

NODAY 

STRANERSTREB 

NEIGHBOURHOOD 

CHECK 

BACKTRACK 

PATH 

CHAINING METHOD 
NEIGHBOURHOOD 
CHECK 
BACKTRACK 
PATH 

INITIi\I.IZATION 
PODINE 
MIDORDERING 
AREA ORDERING 
NEAREST 


Fig, 6,1: Nesting of Procedures 
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Procedure » STEMERSTHBE « (aM ’CHAimS METHOD *>. This 
has been briefly outlined in Sec, 6,4.1 step 9 and 12. 
Therefore this paragraph only reviews the procedures local 
to these two procedures which are also similar. The 
procedure •Neighbourhood' begins with the 'start cell' stored 
in the variable 'LISTl', and fills up the list with all the 
available cells for laying path between it and the end cell 
(or target cell) by means of procedure 'Check'. Procedure 
check passes a cell for inclusion in the list if its 
'path code' is free (i.e. =0) or has the same as that of 
'end cell'. The second condition is important if connection 
is on second side. All four neighbouring cells of each 
new cell added in the list is checked in this manner. Then 
by means of two list pointers, function 'path' checks 
whether a path exists. This returns True only if some cell 
added in the list has the same path code as the path number 
under consideration in the case of 'Straner's tree' and 
if the ce3.1 number is same as that of end cell in the case 
'chaining method*. i.lso each cell added in the list has its 
'link cell’ field updated except the start cell for which 
this field is zero,. If a path exists the procedure 'Back 
track* traverses the cells in the list by means of 'link 
cell’ field (i.e. 'chain coord' field in program) starting 
from the end cell or target cell, it outputs the cell 
numbers when-ever a change of direction is involved. 
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PROCEDURE ' EEAKSST * i Phis takes ■unconnected 
( 'xmtried’) pins one hy one in the selected ro-wr 'of 'net 
matrix* and calculates the nearest already connected pin, 
whose numbers are stored in the 'pin n'umber array*. The 
distance between the two and numbers of the untried pin 
and already tried pin are stored in its local array 'B'. 
Initially the first pin of the net is considered as’tried* 
or already connected. The local array B will confeist 
of as many elements as there are unconnected or untried 
pins in the net. This number is returned in the parameter 
'KJ ’ , The elements of 'B * are then sorted and that 
untried pin which corresponds to the element giving the 
least distance is treated as next 'start cell', After 
invoking the routing procedures this pin is added to the 
'pin number array'. 

6.5 'MIC FILE INSTRUCT I OHS FOR P.R0URAM EXECUTION 

The reasons for MIC file instructions and Appendix 
VII were mentioned in Sec. 5# 4. Tho ’RoutAT* instructions 
for example can be used to execute the router program for 
circuit after the execution of placement program or when 
a placement program output exists in the name of 'PLA1.0UT'. 

If independent input is given, the MIC file instruction can 
be still used by giving it the name 'PLAI (or P1A2 etc). out'. 
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The router program may he followed by 'plotl ' 
for example and then *name1’. Otherwise 'namel ’ inetructions 
must be executed by the simple command *Do Uamel*, 

In the case of plotting routine again renaming 
of files and window and view port dimensions change in 
program data are required. If the necessary input files 
exist this can be carried out independently and then 
'Namel ' can be done. Otherwise this must follow a router 
program. 

See Appendix VII for an example of relevant ‘MIC’ 


file instructions 



CHAPTER 7 


CONCLUSION 

All authors [l»2] agree that there must he a great 
deal more experimental and theoretical investigations 
of the layout techniques. 

7.1 PLACEMENT 

The relative value of each of the method seen 
depends upon the particular application and the amount of 
computer time/storage available. In addition time factor 
limits the variety of test cases that can be employed. 

One good combination would be a fast initial 
placement method seen, coupled with an iterative placement 
improvement method. On the other hand the heuristic 
technique of applying vector forces for placement has the 
advantages: 

(a) Tends to be independent of initial data^ 

(b) Connectability requirements are met along with 

placement . 

(c) The modelling of circuit board is very flexible 
and variety of wiring rules and requirements can 
be satisfied. 

(d) Computational efficiency is satisfactory 

(e) Provides a practical approach and in the cases 
tested does reduce line length considerably. 



91 


7.2 ROUTIIG 

A ntimber of techniques, both, formal and heuristic 
have been developed. Their applicability and usefulness 
vary considerably with the problem under consideration. 

A rough measure of the difficulty of the problem 
is as follows; Suppose ’L' is the total ('Manhattan') 
distance warranted by the interconnections once the wire 
list is determined, and let 'S' be the total number of grid 
segments (cells) available. Then, the ratio S/L is a 
measure of difficulty of the problem. 

like the implementation shown in Chapter 6, most 
are based on lee's Algorithm, The variety stems from 
speedup techniques, net ordering and ordering a pins 
in the net. 

Other techniques that have been tried include a 
graph-to-layout correspondence [9] in which a graph 
structure defines the connections and the router is capable 
of decoding this. Another method [8] assigns all horizontal- 
wires to one side and vertical wires on the other side, 
in the ’channel’ between modulesj^, segments of them are then 
allotted to nets by the router. This requires much less 
computer time and storage but limited in applications. 

It can be said that no routing program can guarantee 
to layout a l 1 interconnections on board (sicde the ratio, 
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S/L may be too small). Hence the possibility of ’cant 
connects* must be allowed for, which must be kept at lowest 
level. 

7.3 SU&GESTIOHS FOR FDTUEB IMPIEMENTATIOH 

Placement ; 

(a) More experimentation to gain use of some stati- 
stical studies in the present implementation. 

(b) Implementation of a monitor that would dynamically 
adjust program parameters such as stopping 
condition etc. depending on circuit, 

(c) Implementation of ’Expand' process so that place- 
ment can be made without regard to predetermined 
sites, 

(d)Use of Pythagorean vectors in sort-delete 
technique of calculating resultant instead of 
X and Y resultant separately as in the present 
implementation. 

(e) Calculation of mobility factor for each module 

instead of single mobility factor for all modules 
as in the present case. 
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Routing ; 

(a) Speed up techniques for Lee’s algorithm. 

(b) Modification to router (and corresponding 
modification to placement to output wire list 
i,e, , cell pairs) to implement all wiring 
rules considered in placement. 

(c) Channel assignment method [8], 

(d) Addition of a software interactive facility to the 
plotting routine for Routing the unmade connections 
(after execution). This can use the cursor 
movement facilities presently available (for 

the Tektronix graphic terminal) with a kind 
of editor added to it. 
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2201 
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03100 

03200 
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03500 
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t 

DIMENSION V{4),*J(4) 

Integer A»8,)(SPAN,7SPAM,XA,irA,xti,in,j,j 

DATA V/0, 0.0, 75. 0.0. 1.9/ * \ 

DATA W/0, 0.60, 0,0. 0,89, 9/ 

ACCEPT *,9 
CALL Nirnevta) 

CALL^ «ilND/<Cw) 

CALL VPOPKO 
CALL SELDEVOJ 
READ(20,»JXSPAM,YSPAN 
PEAD(20,*}A 
IF(A,E0.99998JGJ TO 2 

Y s FLOAT(A)/FLOaT (XSPAO 
I=Y 

Xs: A - I*XSPAN 

V = YSPAM -I 

call LI.^fE (X,Y+0,2,0) 

CALL LTNF (X,Y-0.2.1) 

GO TO J 
ReAUC21,*JA 

IFrA.E0.99999,QH.A.E0.99997jGJ TO 3 

IK (A.FQ,99996) GO TO 2 

READ{2l,*JR 

IK (B,EQ.O)Gn TO 2 

YA = FLtJA'tCAJ/FLOATlXSPAN) 

XA = A - YA+XSPAN 

YA = YSPA,^ - YA 

Yb aFL0AT(8)/FL0AT(XSPA?O 

XB a B - YB^XSPAN 

Yb a ySPA?J-XB 

WRITE(24,4),XA.YA.XB,YB 

CALL LTNI CXA.YA.O) 

CALL LlNItXb,YB,l) 

AsB 

GO TO 4 
CONTINUE 
ACCEPT*, J 

IFfA, £0.99997) GO TO 2 . 

STOP 

END 



SQ PLACE. PAS 
R 60 

XSPAN * 60f YSPAN aBO?. 

R 950 

*lC2l :* 60,0? VC2J :» 0,75; 

R 970 

«t4l s» 80 . 0 ; vt4] :« 0,75? 

£d 

REMAME INPOTsPLAl.iN 
R SPAS 
PLACE. PAS 
ASS TTY 3 

ex PLACEfSYS:PASLNK./5EA SYSSGPGi 


^FORLiIB 



RDuuoui 
Dfili Dun 
Dgfj PUAl.OUT 
RENAME PbAt.lNsINPUT 
RENAME IMPUr=OtJTPUT 
R PASREb 
ROUrE.PAS 
EX ROUTE .REb 

RENAME PLAI.OUTsINPUT 
RENAME ROUl.aUTsOUTPlIT 
RENAME. OUTlsOUT 



00100 

00200 

00300 

00400 

00500 

00600 

00700 

OOftoO 

00900 

01000 


SU PLOT, FOP 
R 300 

* daTA V/0,0, 0.75, 0.0, 1.0/ 

R 4o0 

♦ daTA W/0. 0,60, 0,0. 0,80.0/ 
EB 

rename FnR20.UAT=PTjAl.OUT 
RENAME FOH21 .OAl=8aul .OUT 
ASS TTV 3 

EX PLUT.FOR,/SFA SYSSGPGS 



90100 

00200 

00300 

00400 


rename pJal ,out=f or20,dat 

opQflS? >;?yl«$utafor21 .dat 
REVAMfc Pl.Al ,TN = lMp'ir 
DjSL OUTpUl 



Exit 
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